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We consider a process to create quasi long-range quantum discord between the non-interacting end spins of a 
quantum spin chain, with the end spins weakly coupled to the bulk of the chain. The process is not only capable 
of creating long-range quantum correlation but the latter remains frozen, when certain weak end-couplings 
are adiabatically varied below certain thresholds. We term this phenomenon as adiabatic freezing of quantum 
correlation. We observe that the freezing is robust to moderate thermal fluctuations and is intrinsically related 
to the cooperative properties of the quantum spin chain. In particular, we find that the energy gap of the system 
remains frozen for these adiabatic variations, and moreover, considering the end spins as probes, we show that 
the interval of freezing can detect the anisotropy transition in quantum XY spin chains. Importantly, the adiabatic 
freezing of long-range quantum correlations can be simulated with contemporary experimental techniques. 


I. INTRODUCTION 

Quantum correlation [1, 2] is one of the principal charac¬ 
teristics that separates the quantum domain from its classical 
counterpart. The properties and phenomena that arise from it 
are neither reproducible nor simulatable in classical systems. 
The amount of quantum correlation that exists between two 
subsystems of a pure quantum state is completely captured 
by entanglement [1]. However, for mixed states, local mea¬ 
surements may reveal nonclassical features, that are present 
even in non-entangled or separable states [2]. Measures that 
capture quantum correlations beyond entanglement, such as 
quantum discord (QD) [3], have been instrumental in inves¬ 
tigating several protocols of quantum information and com¬ 
putation (QIC) [4-6], quantum phase transitions [7], many- 
body dynamics [8, 9], quantum biology [10], and metrology 
[11, 12] (see [2], for a review). 

Methods to create long-range entanglement have attracted a 
lot of attention due to their importance in several protocols in 
QIC [13-20]. Such investigations led to the discovery of pro¬ 
cesses like entanglement swapping [13] and repeaters [14], 
and concepts like localizable entanglement [15]. In the last 
decade, there have been several instances where quantum cor¬ 
relation measures such as QD were claimed to be important 
[2, 4-12]. For example, it provides an interesting perspective 
on the dynamics of open quantum systems [21], where entan¬ 
glement is fragile. In quantum systems subject to noisy envi¬ 
ronments, for both Markovian and non-Markovian evolution, 
QD is more robust than entanglement [22, 23]. Interestingly, 
for certain types of quantum states under dephasing, QD ex¬ 
hibits a qualitatively different robustness. It remains/rozen for 
finite evolution times [24], even though entanglement suffers 
sudden death. In recent times, a significant amount of research 
has been undertaken to characterize the phenomena of freez¬ 
ing of quantum correlations [25-27], including entanglement 
[28]. 

In this letter, we report quasi long-range QD between the 
end spins of a finite quantum spin chain, when the end spins 
are weakly coupled to the bulk of the spins. The work is mo¬ 
tivated by the desire to investigate the quantum correlation 
properties of the end spins due to the collective effect of the 
bulk spins that act as a reservoir for the end spins, thus giving 
rise to interesting long-range quantum phenomena. In partic¬ 


ular, we find that the long-range QD of the ground state can 
exhibit non-temporal freezing, while the weak end-couplings 
are varied. We term the freezing as “adiabatic”, as the phe¬ 
nomenon is observed while the weak end-couplings are adia¬ 
batically varied. No such freezing behavior is observed for the 
long-range entanglement. The observed phenomena makes 
long-range QD a robust resource, against moderate thermal 
fluctuations and variations in certain system parameters, for 
implementing quantum protocols between distant spin qubits. 

Interestingly, and in contradistinction to temporal freezing 
[24-27], the observed adiabatic phenomenon is characteristic 
of the quantum system rather than an external environment, 
and is intrinsically related to the many-body properties of the 
quantum spin chain. For instance, in an experimental setting, 
one may consider the two end-spins to be probe sites weakly 
coupled to a system consisting of a spin chain [17-19], where 
the weak end-couplings can be controlled. The long-range 
QD between the probe spins can then be shown to identify the 
nature of interaction in the system. Specifically, the freezing 
interval of long-range QD can be used to define an order pa¬ 
rameter that detects the “anisotropy transition” in the system 
considered. Further, one can show that the adiabatic freezing 
phenomena is not limited to long-range quantum correlations 
but is also manifested in other cooperative properties of the 
quantum spin chain, such as the energy gap. 

We note that there are recent experimental proposals and 
techniques to generate and characterize long-range quan¬ 
tum correlation in quantum spin systems [29-31] (also see 
Refs. [32]). In this light, we investigate the effect of ther¬ 
mal fluctuations in the spin system, which is important for 
potential experimental implementations, and And that the adi¬ 
abatic freezing is stable below a critical temperature. Hence, 
our results show that in an experimentally accessible regime 
of weak end-couplings and low-temperature, high quasi-long 
range quantum correlations can be adiabatically frozen for ap¬ 
plications in various quantum information protocols. 


II. METHODOLOGY 

Let us consider an anistropic XY quantum spin chain con¬ 
taining N spins with a closed end. The Hamiltonian for such 
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a system can be written as 

N 

i 

where cr^^i = Ji and Ki are the dimensionless inter¬ 

action strengths. k(> 0) has the unit of energy, a'^, i = x,y, z, 
are the Pauli spin matrices. The open-end case is obtained by 
setting Jpf and ICn equal to zero. 

The quantum correlation between two arbitrary sites, in 
the ground and thermal equilibrium state of the Hamiltonian 
can be obtained by deriving the two-site reduced density ma¬ 
trix, following the seminal work in [33]. The symmetry of 
the Hamiltonian ensures that all (single-site) magnetizations 
(ct“), V a = (x, y, z) vanish in the absence of any external 
fields. The only non-vanishing two-site terms are (crfcrj), 
(afcrj), and (crftT|). Hence, any two-site reduced density 
matrix, for arbitrary sites i and j, can be written as pij = 
1/4(1 + J2c=x,y,z are 

the two-site correlation functions and I is the two-qubit iden¬ 
tity matrix. Tj““ can be analytically derived by solving the 
Hamiltonian in Eq. (1) [33, 34]. 

Since pij is Bell-diagonal, its QD [3] can be analytically 
obtained [35]. The eigenvalues of pij, in terms of Tj““, is 
given by {ej equal to 1/4(1 ± and 

1/4(1 ± - Tf^) + T,f), for i = 1 to 4. The quan¬ 

tum mutual information is given by the relation, I{pij) = 
Ci log2(4ei), and the classical correlation, obtained after 
optimization over measurements on a single-party, is given by 
C{Pij) = J2l=i Xk\og 2 { 2 xk), where Xk = {I + (-1)^ x)/2 
(fc = 1,2), and X = max{|7;7|, \T^/\, \^^\}. The QD is 
then given by the relation, 'D{pij) = I{pij) — C{pij) = 

Yh=i log2(4ei) - LLi log2(2xfe). Similarly, the en¬ 
tanglement, using concurrence [36], between any two sites 
can be analytically derived in terms of the correlation func¬ 
tion [34]: £{pij) = max [0, 2 max[{ei}] — 1]. 

III. ADIABATIC FREEZING OE QUANTUM DISCORD 

Let us consider an A^-spin non-periodic quantum spin 
chain, with nearest neighbor interactions, where the two spins 
at the edge of the chain (end spins) are weakly coupled to the 
remaining bulk of — 2 spins. The Hamiltonian for the spin 
is given by 'H = (k/4) {Ubuik + 'Hend), where 

-Hduik = «+i + afaf+i), (2) 

Uend = Ai(erferf -f + A 2 (cr/erf -I- er](,_^er](,)(3) 

The above Hamiltonian can be exactly solved to calculate the 
quantum correlation between any two spins in the chain. The 
cases for which Ai = A 2 < Ji, ICi, Vi = 2,..., N — 2, will be 
referred to as “balanced” weak-coupling. It is known that for 
Ai(= A 2 ) —s- 0, the two end spins in the ground state of the 
chain are maximally entangled. Denoting concurrence [36] 
by £l, we have £l I'm that limit. This gives rise to quasi 
long-range entanglement, interestingly, in a quantum system 


with very short-range interaction [16, 17]. In the balanced 
case, decreases monotonically with increasing Ai(= A 2 ), 
and the rate of decay increases with system size. Also, £l 
may exhibit non-temporal sudden death (cf. [9]). 



FIG. 1. (Color online.) Adiabatic freezing of quasi long-range quan¬ 
tum discord ('Dl)- The bulk is an XX spin chain. For each non-zero 
value of A 2 , freezing of TPl takes place for adiabatic change of Ai. 
The behavior of entanglement is given in the inset. All quantities 
are dimensionless, except QD (in bits) and entanglement (in ebits). 
Here, N = 10. 


Let us now consider the case where T-Lbuik represents an 
ordered, XX spin chain, such that Ji = ICi = 1,V i = 
2,N — 2, and the coupling strengths at the ends are such 
that Ai, A 2 < Ji- When Ai J A 2 , we refer to it as “unbal¬ 
anced” weak-coupling. Under such coupling, the behavior of 
and long-range QD {T>i) i.e., the QD of the reduced den¬ 
sity matrix of the end spins, piN, of the ground state, are qual¬ 
itatively different. Specifically, for fixed values of A 2 and on 
slowly increasing the parameter Ai, from approximately 0 to 
1, Dl freezes, i.e., remains unchanged in value for the range 0 
< Ai < A 2 . The value of the frozen Dl (^l) dependent on 
the fixed A 2 and the size of the spin chain, N. However, the 
freezing interval {I j), i.e., the region on the Ai-axis over which 
Dl remains frozen, is equal to A 2 . As Ai is increased beyond 
A 2 , i.e., for the range Ai > A 2 , the freezing of Dl ceases and 
it decays exactly similar to the balanced case. Interestingly, 
the behavior of entanglement does not distinguish the regimes 
0 < Ai < A 2 and Ai > A 2 . In particular, £l decays in a 
similar fashion in both the regimes, though the maximum £l 
at Ai ~ 0 and the value of Ai at the non-temporal death (i.e., 
Af*, such that for all Ai > Af, £l = 0), both decrease with 
increasing A 2 . Ligure 1 shows the freezing of long-range QD 
in an A^ = 10 spin XX chain for different values of the fixed 
weak coupling, A 2 , with the variation of Ai. We call such phe¬ 
nomenon as “adiabatic freezing” since it can be observed in 
the adiabatic evolution obtained by slow variation of the weak 
coupling Ai. The phenomenon of adiabatic freezing can be 
explained by observing the behavior of the two-site correla¬ 
tion functions that are obtained by solving the Hamiltonian in 
Eq. (3). Note that the freezing of Dl observed in the present 
context is purely a property of the quantum spin chain, and 
not a conjunction of the system and an environment. 

We note that both and Dl are quasi long-range [17], i.e.. 
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their value decreases with increase in system-size N. Hence, 
for a fixed A 2 , the frozen discord value, V^, decreases as N 
grows larger. However, fixing the weak-end couplings to low 
values, say Xij Ji ~ 0.001, for N=2QQ, a reasonably high 
Ri 0.93 can be obtained, and frozen upto one order higher, 
i.e.. If = \2 = 0.01. Moreover, even symmetric quantum cor¬ 
relation measures, such as symmetric QD [37], exhibits adia¬ 
batic freezing for Ai 7 ^ A 2 . 





FIG. 2. (Color online.) Variation of X(piAr), C{p\n), and with 
increase in the end coupling Ai, for a fixed A 2 = 0.2, and for N = 20 
spins, where the bulk is an XX chain. T>l is frozen for Ai < A 2 . 

= 0.322 and I / = 0.2. All quantities used are dimensionless except 
Vl,I{pin), and C(pijv), which are in bits. The inset magnifies the 
red-encircled region. 


The analysis of two-site correlation functions between the 
end spins, can shed light on the observed adiabatic freez¬ 
ing of QD. From Fig. 2, for N = 20 spins with A 2 = 0.2, we 
observe that < 0, and \T-^^\ < 1, V a. Since, = 
- T™ [33, 34], therefore and are the only in¬ 
dependent variables, with \T^^\ < \T^^ \ and \Tl^\ < \T^^\. 
Moreover, it is seen that is constant during the adiabatic 
evolution. Hence, the quantum mutual information (I(pijv)) 
is an entropic function of 1 ± Since, |rf^| > \Tf^\ > 
\T^^\ for Ai < A 2 , the classical correlation (C(piAr)) is also a 
function of 1 ± Hence, I(piAr) and C(piAr) decay with 
identical rates leading to a constant long-range QD. However, 
for Ai > A 2 , \T^m\ > \I'in\ > \TIn\ C{pin) is now a 
function of 1 ± Since, is constant during the evo¬ 
lution, C{pim) freezes for Ai > A 2 , in contrast to I(pijv) in 
that range. Hence, the freezing of the long-range T>l{pim) 
disappears for Ai > A 2 . As A 2 —{0, Ai}, there is no freez¬ 
ing of Di(pi 7 v)j as C(pi 7 v) is always constant. In contrast, 
for entanglement, no adiabatic freezing occurs. In Fig. 2, 
since 0 < \T^^\ < 1, V a, £1^ is given by the relation, £l 
= max [0, l/2{giN - Hin)], where, giN = |Tf^ + T^^\, and 
hiN = 1 — \Ti§\\T^^\ > 0. As Ai increases gi^ and hiN 
decreases, due to decreasing \T^^\, and constant \T^^\. For 
giN > hiM, £l decreases with Ai, and at giM < hiM, £l = 
0 (sudden death). For Ai = A 2 , we have = Tfff = z (say). 
Therefore, Sl = max [O, l/2(z^ -F 2\z\ — 1))] [34]. 

The observed behavior of I{pi]^) and C{pi]\f) is consistent 
to what is observed in the freezing of QD, first observed in 
dephasing quantum systems [24], that heralded a substantial 


amount of research in recent years [25-27]. However, we note 
that the freezing of 'Dl{pim) observed in the present context 
is purely a property of the long-range correlations in the spin 
chain and is devoid of any external decoherence, in contrast to 
conventional freezing phenomena. 


IV. DETECTING ANISOTROPY TRANSITION 

The adiabatic freezing of the quasi long-range QD (Vl) in 
spin chains, with Hamiltonians of the form given by Eq. (3), 
can be used to investigate certain intrinsic properties of the 
bulk Hamiltonian, such as the anisotropy transition. Though 
there exist methods to detect anisotropy in quantum spin sys¬ 
tems [38], our results provide an information-theoretic per¬ 
spective to investigate these quantum properties in an experi¬ 
mentally viable manner. 

Let T-Lbuik be an ordered XY spin chain with Ji = 1 + 7 
and /Ci = 1 — 7 , V z = 2,..., N — 2, with 7 ( 7 ! 0) being the 
finite anisotropy present in T-Lbuik- The end-coupling is unbal¬ 
anced, such that Ai 7 ^ A 2 . We focus on the parameter regime 
Ai, A 2 < 1, where the relevant physics under consideration 
is clearly observed. Let us now consider the adiabatic freez¬ 
ing of Vl in the above spin chain, with possible anisotropy in 
T~Lbuik, by varying the weak end-couplings. The presence of 
anisotropy introduces changes in the freezing characteristics 
of 27^. Lor instance, increases with 7 upto a certain fixed 
value that depends on N, and then decreases. See inset of 
Lig. 3. Interestingly, the freezing interval, I f undergoes a tran- 



FIG. 3. (Color online.) Detecting the anisotropy ( 7 ) transition. The 
order parameter, A 2 — (/ is plotted against 7 , with fixed A 2 . The inset 
shows that can be substantially increased for positive 7 , below a 
threshold value. The anisotropy transition can also be seen from the 
change in the curvature (convexity to concavity) of with respect 
to 7 , where the A-dependent critical values approach 7 = 0 through 
negative 7 . All quantities are dimensionless, except QD (in bits). 


sition that can detect the anisotropy transition at 7 = 0 in the 
bulk spin system. This allows us to define an anisotropy tran¬ 
sition order in terms of the system parameters and the freezing 
interval If. This can be seen as follows. We evaluate If for 
different values of Ai, and a fixed A 2 . By using A 2 — Q as the 
order parameter, one can detect the 7 -transition in the sys¬ 
tem at 7 = 0 , driven by changes in the anisotropy parameter. 
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The order parameter X 2 — If is finite for anisotropic systems 
and vanishes at 7 = 0. Further, one observes that X 2 — If 
sharply decreases as 7 is increased beyond an A^-dependent 
threshold value (7^). It is observed that as N is increased 
upto certain values, the critical parameter 7 ^ —0. Hence, 
for large but finite N, transition in A 2 — // at captures the 
7 -transition of the class of quantum spin models described by 
the Hamiltonian in Eq. (3). The observed transition is shown 
and described in Fig. 3. 

In an experimental setting, one may consider the end spins 
to be probe sites at the edge of a quantum spin chain con¬ 
sisting of the bulk spins. The two end-spins can be defects 
or scattering particles in a spin chain [17, 39], and the weak 
end-coupling can be experimentally controlled. Under such 
conditions, the anisotropy of Hbuik can be detected using the 
order parameter X2 — If, where the desired quantities are ob¬ 
served from the freezing of long-range QD between the two 
probe spins. Hence, the quantity A 2 — Z/ can serve as a suitable 
order parameter in many-body simulations. The advantage of 
the approach lies in the fact that all performed measurements 
and tuning of interactions are associated only with the probe 
spins that weakly interact with the quantum spin chain. The 
detection of the intrinsic anisotropy in the bulk spin chain is 
achieved without disturbing the system. This provides an in¬ 
teresting role for adiabatic freezing of discord in investigating 
many-body phenomena. 

V. FREEZING OE ENERGY GAP 

Another interesting phenomenon that connects adiabatic 
freezing of long-range quantum correlations to the coopera¬ 
tive properties of the quantum spin chain, is the freezing of en¬ 
ergy gap. For the considered model, the weak end-couplings 
between the bulk and end spins introduces a finite energy gap. 



FIG. 4. (Color online.) Adiabatic freezing of the energy gap. Nu¬ 
merically obtained values of Ag is plotted against Ai, for different 
fixed values of A 2 . The spin model is the same as in Fig. 1, but 
for N = 20 (red) and 30 (black). For Ai > A 2 , Ag is frozen. The 
difference in frozen Ag for different N, marked in one case with a 
double-headed arrow, increases with A 2 . All quantities plotted are 
dimensionless. Note that the behavior of Ag and T>l with response 
to Ai, for fixed A 2 , is complementary. 


The energy gap is an intrinsic property of a quantum spin 


system that can be obtained from the excitation energy spec¬ 
trum or the dispersion relation of the function [33, 34]. 
Let us consider the system, where T-Lbuik in Eq. (2), represents 
an ordered XX spin chain. The only anisotropy in the system 
arises from the coupling strengths at the end spins. Solving 
the Hamiltonian, one can obtain the excitation spectrum, A^. 

For the balanced case, Ai = A 2 = A, the dispersion relation 
of the excitation energy is given by Afc = cos(A:), where k are 
the quasimomenta modes. These modes satisfy the eigenvalue 
equation [16], pcot(A:)[cot((A — l)fc/2)]'^ = A^/(2 — A^), 
for A 7 ^ 1, with jj, = ±1 being the eigenstate parity. Consid¬ 
ering /r = 1, the energy gap is then given by k', which min¬ 
imizes Afc/ = cos(fc'), while satisfying the above eigenvalue 
equation. For A 3> 1, at fc' = tt/2 — 5, where (5 —0, an an¬ 
alytical expression for the energy gap (Ag) is obtained [17], 
where Ag ( 7 r/ 2 A)[(l + 2 /(A(AV (2 - A^) + 2 )] = /(A). 
Numerical analysis shows that for the case, Ai 7 ^ A 2 , the 
quasimomenta k' corresponding to the energy gap satisfies the 
eigenvalue equation for A = min[Ai, A 2 ]. Therefore, the en¬ 
ergy gap is given by the relation, Ag = min[/(Ai), /(A 2 )], 
where /(A), defined above, is a monotonically decreasing 
function of A. The adiabatic freezing of the numerically es¬ 
timated energy gap, is shown in Fig. 4. Using the above 
relation for Ag and monotonicity of the function /(A), it is 
obvious that Ag freezes for A 2 < Ai, since in this region, 
/(A 2 ) < f{M) and Ag is thus independent of Ai. 



FIG. 5. (Color online.) Variation of the analytically and nu¬ 
merically estimated values of the energy gap with increasing size 
of the spin chain. The analytical expression is given hy Ag = 
min[/(Ai),/(A 2 )] (blue-square), and is observed to be consistent 
with numerically obtained values of Ag (red-circle) at large N, for 
Ai = 0.4 and A 2 = 0.6 (and, Ai = 0.6 and A 2 = 0.4). The inset shows 
T>l (maroon-diamonds) with increasing N. The figure shows that at 
large N, Ag and 'Dl scale as 1/N. 

The adiabatic freezing of Ag is complementary to that of 
in terms of the variation of weak end-couplings. Figs. 1 
and 4 show that while Vl freezes for Ai < A 2 , Ag is constant 
for Ai > A 2 . Moreover, while the frozen value of Vl de¬ 
creases with increase in A 2 , it does the opposite for Ag. Fig¬ 
ure 5 shows the agreement between the analytically and nu¬ 
merically obtained values of Ag, at large N. We observe that 
at large N, the frozen values of both Ag and 17^ scale with 
1 /N [34]. The behavior of Ag shows that the phenomenon of 
adiabatic freezing is manifested through the intrinsic proper- 
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ties of the quantum spin chain. 

VI. RESPONSE TO THERMAL ELUCTUATIONS 

In the previous sections, the phenomena of freezing is ob¬ 
served for the ground state of the quantum spin chain under 
consideration. For small thermal fluctuations, the system is 
no longer in the ground-state, but a mixed state in equilibrium 
at some small temperature (T). To obtain the reduced two- 
site density matrix of the thermal equilibrium state of the spin 
system, at temperature T, one must find the thermal two-site 
correlation functions Tj““(/3) [33, 34]. Here, (3 = l/ksT, 
where ks is, the Boltzmann constant. 



kj/x 

D 


FIG. 6. (Color online.) Response of freezing to thermal fluctuations. 
We plot the ratio of the thermal freezing interval (l^) and the ground 
state freezing interval (If) as a function of temperature. Here, X2= 
0.2. We observe that for every N, there exists a critical temperature 
(T ^) beyond which adiabatic freezing disappears for finite chains. 

decreases linearly with N, as shown in the inset. The solid 
green line, in the inset, is the linear fit, y = 19.72 — 0.184®. All 
quantities plotted are dimensionless.The abscissae in the figure and 
the ordinates in the inset are multiplied by lO'*. 


We observe that adiabatic freezing persists at finite T, 
below a critical temperature T^, that also depends on the 
system-size N. The value of frozen discord is constant 
for temperatures below , though the freezing interval de¬ 
creases as T increases. is a linearly decreasing function of 
N. Figure 6 shows the effect of thermal excitations on freez¬ 
ing interval of long-range QD in an XX spin chain. 

VIE DISCUSSION 

In recent years, experimental generation of long-range 
quantum correlation in antiferromagnetic spins chains have 
been reported using strontium-cuprate compounds, such as 
Sri4Cu2404i, to simulate dimerized spin-chains [29]. Quan¬ 
tum correlations are experimentally measured using low tem¬ 
perature magnetization, magnetic susceptibility, and heat ca¬ 


pacity [29, 32]. Another experimental protocol, relevant to 
the quantum spin chain considered in our study, was proposed 
using ions of the ytterbium isotope, ^^^Yb, in segmented lin¬ 
ear Paul traps [30], utilizing the tailoring of axial trapping 
potential to generate the spin-spin coupling, and using mi¬ 
crowave pulses to generate effective spin interactions. A more 
recent study devises a similar scheme using superconducting 
flux qubits, using dc currents and microwave pulses to control 
the spin-spin interaction [31]. Moreover, recent experimental 
studies have also observed the freezing phenomena in various 
quantum systems [25]. 

In our work, we And that there exists adiabatic freezing of 
quasi long-range quantum correlations in finite quantum spin 
chains. We show that the observed phenomena is robust to 
the weak end-spin couplings and finite thermal fluctuations, 
which are the fundamental elements in experimental control 
of quantum systems. Further, by tuning the end-spin coupling 
one can obtain relatively high values of long-range QD. In¬ 
terestingly, we have observed that a finite interaction between 
the two end spins can vastly increase the shared QD between 
the sites, without affecting the freezing interval. This is intu¬ 
itively plausible, as the finite end-to-end coupling encourages 
greater correlation between the end spins. Alternatively, one 
can also study the phenomena in the XX and XY spin chains 
in a transverse magnetic held. Preliminary investigations re¬ 
veal that, for weak end-couplings, an effective freezing [27] 
of long-range QD can be characterized. For quantum spin 
chains such as the frustrated spin-1/2 Ji — J 2 model, spin-1/2 
XX chain with alternating interactions, and the spin-1 AKLT 
model, long-range quantum correlations are observed [17], al¬ 
though adiabatic freezing is absent. We also note that the phe¬ 
nomenon of adiabatic freezing can be utilized to experimen¬ 
tally detect properties of quantum spin-baths, modelled by an 
interacting bulk spin Hamiltonian, and probed by weakly in¬ 
teracting spins at ends of the bath. 

To conclude, we And the phenomenon of adiabatic freezing 
of quasi long-range QD in the closed dynamics of many-body 
quantum systems. Our work makes a connection between the 
temporal freezing of correlations, observed only in damped 
quantum systems, to the feature of long-range correlations in 
quantum spin chains. However, in contrast to temporal freez¬ 
ing, the adiabatic phenomena is an intrinsic property of the 
considered spin system. It has the ability to detect important 
cooperative phenomena in quantum spin models and in partic¬ 
ular serve as an order parameter for detecting the anisotropy 
transition in quantum XY models. We note that the phenom¬ 
ena is also observed for other system properties such as energy 
gap, and other quantum correlation measures such as symmet¬ 
ric discord [37] and one-way quantum work-deficit [40]. 

With unprecedented developments to simulate quantum 
spin chains in different physical substrates and experimen¬ 
tal techniques to characterize quantum correlations, the phe¬ 
nomenon of adiabatic freezing allows the generation of robust 
long-range quantum correlations between distant parties, for 
application in future quantum technologies. 
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Supplementary Material 


SI. SPIN HAMILTONIAN AND CORRELATION 
FUNCTIONS 


Let us consider an anistropic XY quantum spin chain con¬ 
taining N spins with a closed end. The Hamiltonian for such 
a system can be written as 

N 

■H = ^ ^ (J* «+i -f JCi crfcrf+i), (Sl) 

where cr^^i = Ji and K-i are the dimensionless inter¬ 

action strengths. k{> 0) has the unit of energy, cr®, i = x,y, z, 
are the Pauli spin matrices. The open-end case is obtained by 
setting Jpf and ICn equal to zero. 

The two-site quantum correlation, between arbitrary sites, 
in the ground state of the Hamiltonian can be obtained by de¬ 
riving the two-site reduced density matrix, following the sem¬ 
inal work in [SI]. The Hamiltonian given in Eq. (sl), can be 
transformed in terms of spin-raising and -lowering operators, 
a] = af + iaf and ai = af — iaf, to obtain 


where Jl = {Ji + K.i)l2 and /C- = {Ji - ICi) 12. The partly- 
Fermi, partly-Bose operators (a^) can be transformed to a 
set of strictly Fermi operators (fc^), using the Jordan-Wigner 


transformations [S2], such that ki = exp 




and k] 


a] exp 


Equation (sl) takes 


E l —I ^ T 

j=i a)aj 

a quadratic form in terms of the Fermi creation (fc^) and 
annihilation (k) operators, that can be diagonalized. The 
quadratic-form Hamiltonian is given by H = K^^ii kj-^ijkj 


+ ^{kjBij k] + h.c.), where Aij = ^{J' 6 i+ij + 
is a symmetric matrix and Btj = is an 

anti-symmetric matrix. For closed-ended chains, Ain = Ani 
= J'^ and Bin = —Bni = K,'n- As shown in [Sl], any two- 
site reduced density matrix of the ground state, for arbitrary 
sites, can be derived in terms of the matrices A and B, by 
solving the eigenvalue equation, (j)k{A - B){A + B) = 

The dispersion relation of the function gives us the exci¬ 
tation spectrum that can be used to estimate the energy gap 
in the system. A corresponding vector, ipk, is defined as ipk 
= -^{A + B)(f>k. A unitary correlation matrix, Q is then ob¬ 


tained by the relation, Qij = — J 2 k V'fc* 't’kr 

The two-site reduced density matrix can be obtained from 
the single-site magnetizations, (erf), and the two-site correla¬ 
tion functions, (erfer^) — (erf)(cr^), where (a,/3 = x,y,z). 
The symmetry of the Hamiltonian ensures that the only non¬ 
vanishing (single-site) magnetization is (erf). However, for 
no external fields, the single-site magnetization also vanishes. 
Moreover, the only non-vanishing two-site terms are (erf trj), 
(aferj), and (tTfCT|), which gives us the two-site correlation 
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FIG. sl. (Color online.) The set-up. The end spins are weakly cou¬ 
pled (black-dashed lines) with the bulk, which are strongly coupled 
to each other (black solid line). 


functions, since (erf) = 0, V a. In terms of the correlation ma¬ 
trix, 0, derived in the main text by diagonalizing the Hamilto¬ 
nian in Eq. (sl), the correlation functions are given by [Sl], 



Gi,i+i 

Gi,j 

rpXX _ 



^ij — 

Gj-l,i+l 

■■ 

II 

Gi+l.i 

Gi+i,j-i 

Gj,i 


rpZZ _ 

-^ij — 

{Gi,iGj,j — 

Gi,jGj,i)- 


Note that and are minors of the determinant of Q. For 
nearest neighbors, the above correlations reduce to Tfffi = 
= Qi+i^„ and = -Gi^i+i since 

(erf) = —Qii = 0. For the end-to-end spin (see Fig. sl), two-site 
correlation function, (Tf fr), the minor isanA—IxA — 1 
matrix and the expressions in Eq. (s3) can be simplified to 
obtain the following: 

Ti,% = -Sn.i det {A - B)/\ det(^ - B)\ = Gn.u 

= -Gun det (.4 - B)/\ det(^ - B)\ = Gi,n, and 
T!^n = -Qun Gn,i = T^n- (s4) 

Now, any two-site reduced density matrix, for arbitrary sites 
i and j, can be written as 

= E (s5) 

oc—x,y,z 

where I is the two-qubit identity matrix. Since the derivation 
of the two-site density matrix depends on the diagonalization 
of N X N matrices, such as A and B, the method can be ex¬ 
ecuted for chains with a large number of spins, and allows us 
to study the asymptotic behavior of several system properties. 
In our case, these are the long-range quantum correlations. 

To obtain the reduced two-site density matrix of the thermal 
equilibrium state of the quantum spin system, at temperature 
T, one must find the thermal correlation matrix GijiP), in a 
similar fashion to the analytical derivation done above. Here, 
/3 = l/ksT, where ks is the Boltzmann constant. The corre¬ 
lation matrix can be written as 

0u(/3) = - ^kj{{r]lr]k )0 - {Vkvi)fi), (s6) 

k 
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where ril are the spinless fermionic operators that diagonal¬ 
ize the Hamiltonian, thus generating fermionic excitations 
in the ground states with energy |Afc|. From Fermi statis¬ 
tics, (?7fc»7fc)/3 = l/(exp[/3Afe] -f 1), and hence, Gij{l3) = 
- tanh[/?Afc/2](^fej. Using Giji/3), the reduced two- 

site density matrix for the thermal equilibrium state can then 
be evaluated by following the expressions for the ground state 
of the spin system as shown in Eqs. (s3 - s4). 


S2. MEASURES OE QUANTUM CORRELATION 

The correlation functions Tj““ can be used to derive quan¬ 
tum correlation measures, such as quantum discord and en¬ 
tanglement, for reduced two-site density matrices of both the 
ground and thermal equilibrium states of the system. The ob¬ 
tained two-site reduced density matrix, in Eq. (s5), is Bell- 
diagonal and hence its quantum discord [S3] can be calculated 
using an analytical optimization [S4]. Eor the Bell-diagonal 
density matrix, its eigenvalues, e^, can be obtained in 
terms of 


62 = 1/4(1-7;7 + T™+i;-); 

63 = 1/4(1+ 1/7-7/7 

64 = 1/4(1 + 7/7 + 7/7 - 7/7). (s7) 

The quantum mutual information is given by the relation, 
I{pij) = 6* log2(4ei). The classical correlation obtained 

after optimization over measurements on a single-party, is 
given by the relation 


of the obtained two-site reduced density matrix is given by 


— max [0,2 Cmax 1] 


= max 




,(sl3) 


where emax = max[{677i]’ gfj = ± , and = 

1 ± Tf/. Hence, once the correlation functions, are 

‘‘J •'J 

known from Eq. (s3), quantum discord and entanglement can 
be obtained using Eqs. (sl3) and (si 1), respectively. 

To highlight the role of the correlation functions in the be¬ 
havior of quantum discord and entanglement during the phe¬ 
nomena of adiabatic freezing, we consider an explicit exam¬ 
ple. Let us study the model considered in the main text; an 
A-spin open quantum spin chain, with nearest neighbor inter¬ 
actions, with two spins at the edge of the chain (end spins) are 
weakly coupled to the remaining bulk of N—2 spins, as shown 
in Eig. si. As presented in the main text, the Hamiltonian for 
such a spin chain is given by. 


"H = Hbuik + 'Hend, where 

N-2 

'Hbuik «+i + (sl4) 

i^2 

He„d = J[Ai(a[=af+a7_ia^) 

+ A2(crfcr| + cr7-l'^w)] • 

Ai and A 2 are the weak end-couplings, and {Ji} = {K,i] = 1, 
such that the bulk forms an XX spin chain. Adiabatic freezing 
of long-range quantum correlations is observed, when one of 
the end-couplings (say, A 2 ) is kept fixed, while the other (say, 
Ai) is adiabatically varied. 


2 

C{pij) = ^Xk\og 2 { 2 xk), where (s 8 ) 

k=l 

= (1 + (—l)^a;)/2, forfc = (l,2), and (s9) 
x = max{|7/7|,|T7|,|7/7|}. (slO) 

The quantum discord is then given by the relation, 

4 2 

= ^e*log 2 ( 46 i) - ^a;felog 2 ( 2 xfc). (sll) 

i=l fc=l 

Similarly, using concurrence [S5] as our measure of choice, 
the entanglement between any two sites can be analytically 
derived. Concurrence of a two-qubit density matrix, p^, is 
defined by the relation 

£(py) = max [0, Cl - 62 - 63 - C4] , (sl2) 

where c^’s are the square root of the eigenvalues of the matrix 
Pijplj, arranged in decreasing order, py = erf ® crj p*j erf 0 
erf. Eor the two-site density matrix obtained in Eq. (s5), p/j = 
Pij, and Ci’s are nothing but the eigenvalues of pij, given by 
Eq. (s7), arranged in decreasing order. Hence, the concurrence 



FIG. s2. (Color online.) Variation of long-range two-site correlation 
functions, (green-square), Tf/f (violet-up-triangle), and 
(maroon-down-triangle), with end-coupling strength Ai, for a spin 
chain with N = 20. The other end-coupling, A 2 is kept fixed at 0.2. 
The inset figure shows the behavior of long-range concurrence (blue- 
diamond) and quantum discord (square-red). The adiabatic freezing 
of T>l for Ai < A 2 is evident from the inset figure, whereas en¬ 
tanglement exhibits non-temporal death at Ai = 0.59. All quantities 
used are dimensionless except entanglement (in ebits) and quantum 
discord (in bits). Compare with Fig. 2 in the main text. 

In Fig. s2, one observes the adiabatic freezing of long-range 
quantum discord (Vl) between the end spins, in the ground 
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state of the Hamiltonian defined in Eq. (si5), when the end 
spin coupling satisfies the condition, Ai < A 2 (fixed), V Ai. 
For Ai > A 2 (fixed), Vl decays with increasing Ai. This phe¬ 
nomena is however not observed for long-range entanglement 
{£l)- The behavior of both Vl and Sl upon adiabatically 
varying Ai can be explained through the variation of the cor¬ 
relation functions and shown in Fig. s2. We 

observe that < 0, with \T^^\ < 1 ,\/ a = {x, y, z), for a 
spin chain with N = 20 spins and A 2 fixed at 0.2. Since, = 
jiyy, therefore and are the only independent 
variables, with \T^^\ < |Tf^| and \T^^\ < |Tf^|. Moreover, 
T^§ remains constant with the variation of Ai. 

Consider the region, Ai < A 2 = 0.2, in Fig. s2. We see that 
\Ti^\ > \Ti§\ > Therefore, the C(pij) is dependent 

only on and decreases with increasing Similarly, 

T{pij) is varies with \T^^\ (\T^^\ is constant) and decreases 
with an identical rate, thus allowing Vl to remain frozen for 
Ai < A 2 = 0.2, as observed in the inset of Fig. s2. For 
Ai > A 2 = 0.2, \T-^\ > \Tf^\ > |Tf^|, and C(p,,) is 
dependent only on | \, which is constant. Therefore, C {pij ) 

is constant for Ai > A 2 = 0.2, but the decreasing I{pij) 
forces Vl to decrease, leading to breakdown of freezing. The 
behavior is consistent with that observed in the phenomena of 
temporal freezing [S6]. 

For entanglement, no adiabatic freezing occurs. Since 0 
< l^iTvl ^ V a, the long-range concurrence is given by 
the relation, £l = max [O, l/2{\g'lj^\ — where, 

= |Tf- + Tf^l, and h+^ = I - |Tf-|| > 0. As Ai 

increases ^in decreases, due to the fact that \Tl^\ 

decreases, and is constant. For, £l = 

l/2(|p^^| — decreases with Ai. For 

£l = 0, and long-range entanglement vanishes. The above 
analysis can be compared with known results on quasi long- 
range entanglement. For Ai = A 2 , the system is isotropic, 
and we have = Tf^ = z (say). Moreover, = 
—z^, which gives us the relation for the long-range entan¬ 
glement, £ 1 ^ = max [ 0 , l/2(z^ + 2\z\ — 1))]. Now for x = 
— {SiS^ + SfSlf) = —1/2 2 , the above expression for en¬ 
tanglement reduces to £l = 2, max [O, (x2 -f |x| - 1/4))], as 
shown in [S7]. 


S3. FREEZING OF ENERGY GAP AND QUANTUM 
DISCORD, AT LARGE N 

An important result discussed in the main text of the letter, 
is the adiabatic freezing of the energy gap, and its comple¬ 
mentary relation to the freezing of quantum discord. In par¬ 
ticular, it is shown that for the Hamiltonian in Eq. (sl5), the 
quasi long-range quantum discord between the end spins in 
the ground state of the system is frozen for Ai < A 2 , while 
the energy gap remains constant in the complementary region 

Ai > A 2 . 

Let us briefly describe the presence of energy gap in the 
spin Hamiltonian considered in the study. As mentioned in the 
main text, for the case, Ai = A 2 =A, the dispersion relation of 
the excitation energy is given by Afc = cos(/c), where k is the 
quasimomentum modes. These modes satisfy the following 



FIG. s3. (Color online.) Variation of the analytically and numerically 
estimated values of the energy gap with increasing size of the spin 
chain. The value of energy gap, using the analytical expression, is 
given by Ag = min[/(Ai), /(A2)] (blue-square), where /(A) is de¬ 
fined in Eq. (sl6), and for exact numerical calculations (red-circle), 
for Ai = 0.4 and A2 = 0.6 (and, Ai = 0.6 and A2 = 0.4), for large N. 
The figure shows that at large N, Ag scales linearly as 1/N. 

eigenvalue equation cot(fc)[cot((W — l)fc/2)] = A^/(2 —A^), 
for A 7 ^ 1, where the positive eigenstate parity has been con¬ 
sidered [S7]. The energy gap is then given by k’, which min¬ 
imizes the dispersion relation Afc = cos{k'), while satisfying 
the above eigenvalue equation. 



FIG. s4. (Color online.) Variation of the long-range quantum discord 
(red-circle) with increasing size of the spin chain. The inset shows 
the behavior of the correlation functions, (green-square), 
(maroon-down-triangle), andT/)^ (violet-up-triangle). Ai and A 2 are 
set at 0.01 and 0.1, respectively. The figure shows that at large N, 
TPl scales as 1/A. 

In large N limit, for k' = 7:12 — 5, where b ^ 0, the 
dispersion relation provides an analytical expression for the 
energy gap (A^) as 

^ ^ A(AV(2-A2) + 2) ^ 

Numerical analysis for the case, Ai 7 ^ A 2 , shows that the 
quasimomenta k' corresponding to the energy gap satisfies the 
eigenvalue equation for A = min[Ai, A 2 ]. Thus, the energy gap 
is given by the analytical relation, A^ = min[/(Ai),/(A 2 )], 
where /(A) is defined in Eq. (sl 6 ). Figure s3, shows the 
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agreement between the analytical expression for Ag and ex¬ 
act numerical calculations for large N. One can then show 
that /(Ai) < /(A 2 ), for Ai < A 2 , and /(A 2 ) < /(Ai), for 
A 2 < Ai- The adiabatic freezing of the energy gap, as Ai is 
varied, is thus evident for A 2 < Ai, as Ag is independent of 
Ai in this range. 

It is known that for the model considered in the study, given 


by Eq. (si5), the long-range quantum correlation between the 
end spins is quasi long-range [S7], i.e., the long-range quan¬ 
tum correlation vanishes with increasing N. The behavior of 
Vl with increasing system size is shown in Fig. s4, which 
plots the value of the frozen quantum discord with increasing 
system size. The figure shows that in the large N limit, both 
the frozen energy gap and quasi long-range quantum discord 
scales with 1 /N. 
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